2024.7.23 安定余裕【control】
次のシステムについて例を示す。
$ G(j\omega) = \frac{1}{j\omega + 1}
code:margin1.py
import control as ct
import matplotlib.pyplot as plt
ct.bode_plot(sys1, margins=True)
plt.show()
chk_margin = ct.margin(sys1)
print(chk_margin)
https://scrapbox.io/files/669f517acba0f3001ce793c2.png
margin【control】関数でゲイン余裕、位相余裕、ゲイン交差周波数、位相公差周波数が求められる。
code:result.txt
(inf, inf, nan, nan)
この場合、ゲイン余裕、位相余裕ともに無限大である。